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Abstract 

In this paper we shall study a gauge theory of nodal spinons which appears as a low-energy 
effective theory for antiferromagnetic (AF) Heisenberg models. In most of studies on the nodal 
spinons given so far, the gauge interaction between spinons was assumed weak and nonperturbative 
effects like instantons and vortices were ignored. In the present paper, we shall study strong-coupling 
gauge theory of nodal spinons and reveal its nontrivial phase structure. To this end, we employ 
recently developed lattice gauge theory techniques for studying finite-temperature and finite-density 
gauge theory. At low temperature and low spinon-density region, an AF long-range order exists. As 
temperature and/or density of spinons are increases, a phase transition to nonmagnetic phase takes 
place. Order of the phase transition is of second (first) order for low (high) density region of spinons. 
At a quantum critical point at vanishing temperature T — 0, abrupt change of spinon density occurs 
as a function of the chemical potential. Implications of the results to the heavy-fermion materials 
and the high-T c cuprates are discussed. 
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1 Introduction 



Since the discovery of the high-T c cuprates, search for nontrivial quantum spin liquid for antifer- 
romagnetic (AF) materials is one of the most challenging topics in the condensed matter physics. 
Among them, possibility of appearance of spin ^ fermionic or bosonic quasi-particles called spinons 
has been studied from various point of view. In certain class of the AF Heisenberg models in two 
dimensions (2D), by using slave-particle representation for the quantum spin operators and mean- 
field theory (MFT) like approximation of decoupling fields, one can simple show that spinons with a 
relativistic dispersion appear as low-energy excitations [1]. If the fluctuations of the decoupling field 
are ignored, spin-spin correlation functions exhibit a power-law decay. This type of spin states were 
dubbed algebraic spin liquids (ASL)[2, 3]. 

In the ASL, besides spinons, there exist gapless gauge-field excitations, which are phase degrees 
of freedom of the decoupling field and its fluctuations may destabilize the MF picture of the ASL. 
It is known that the gauge-field effects in 2D is very strong and low-energy excitations are often 
drastically influenced by the gauge interactions. Generally speaking there are two phases in the 
gauge theory; one is confinement phase in which the gauge field fluctuates very strong and charged 
particles by themselves cannot exist as low-energy excitations. The other is deconfincment phase 
(it is furthermore classified into Coulomb and Higgs phases) in which nonperturbative effects of 
the gauge-field fluctuations can be neglected and charged particles and gauge boson itself are low- 
energy excitations. In the discussion of the ASL, the deconfincment phase of the U(l) gauge field is 
assumed. Though there is no Maxwell term (or no plaquette term) of the gauge field in the original 
action, this assumption is justified if the number of gapless spinons N s is very large because of the 
shielding effects of the gapless matters. 1 . However our recent study on the deconfincment phase 
transition by a large number of gapless matter fields indicates that the critical value of N s is rather 
large. In Ref.[7], we numerically investigated multi-flavor CP 1 models on a lattice and found that 
the critical number of two-component complex bosons for the deconfincment transition is 14 which 
is larger than the number of the two-component Dirac fermions appearing in the typical ASL's. If 
the gauge theory of nodal spinons is in a confinement phase, weakly interacting spinons may be a 
good picture in the distance scale of the original lattice spacing but the groundstate and picture 
of the quasi-excitations at long distances are drastically changed by the violent fluctuations of the 
gauge field. 

In this paper, we shall study the strong-coupling gauge theory of nodal spinons from the above 
point of view. We shall clarify its phase structure and quasi-excitations. To this end, we employ 

1 For discussions on the problem of deconfinement phenomenon in gauge theories coupled with gapless matter fields, 
see Refs.[4, 5, 6] 
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recently developed lattice gauge-theory techniques for studying finite-temperature and finite-density 
gauge theory[8]. 

This paper is organized as follows. In Sec. 2, we briefly review the ASL in which relativistic 
Dirac fermions appear. In Sec. 3, we introduce a lattice model describing strong-coupling gauge 
theory of spinons. Symmetries of the model is also discussed. In Sec. 4, an effective action is derived 
by integrating over the gauge field and spinons. In Sec. 5, by using a MFT like approximation we 
clarified the phase structure of the system. Section 6 is devoted for discussion and conclusion. 

2 Algebraic Spin Liquid 

In this section we shall briefly review the ASL. Let us consider S = 1/2 Hcisenberg model on the 
square lattice, 

Haf = J*^2S X ■ S x+i + • • • , (2.1) 

x,i 

where S x is the spin operator at site x — {x\,X2), i is the direction index (sometimes it also denotes 
the unit vector in the «-th direction), and J > is the nearest-neighbor (NN) exchange coupling. 
In the Hamiltonian (2.1), there may be other interaction terms between spins like the next-NN 
exchanges, and the ring exchanges, etc which are represented by the ellipsis. The spin operator S x 
can be represented by the two-component fermionic spinon operator f x = (/ x ,i,/x,2) = {fx,h fx,l) 
as, 

Sx = \&fx, (2.2) 

where a is the Pauli spin matrices and f x must satisfy the local constraint Y^ a =i fxafxa — 1 because 
the magnitude of the quantum spin is 1/2. By substituting Eq.(2.2) into Eq.(2.1) and decoupling 
the quartic terms of f x in Haf by introducing decoupling fields D xi , the decoupled Hamiltonian is 
obtained as follows, 

H dc = -J2 Dxifl +i fx + H.c. + \D xi \ 2 / J +■■■. (2.3) 

xi 

From Eq.(2.3), it is obvious that the phase of D xi , iA xi = \og(D xi /\D xi \), behaves like a U(l) gauge 
field under a gauge transformation, 

fx^e ia *f x , D xl ^e^D xl e- la *. (2.4) 

In the MF approximation, the fluctuations of D xi are totally ignored and D xi is replaced with its 
expectation value D° xi = (D xi ) in Eq.(2.3). 

Low-energy excitations are determined by the pattern and symmetry of the MF D xi . For example 
in the 7r-flux state, phase of the product of four D x is around the smallest plaqucttes is it. In most of 
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general flux states, the low-energy excitations are described by two-component (or four-component) 
relativistic Dirac fields ip a where a — 1, • • • ,nj is the flavor index (the number of nodes). For the 
7r-flux state, -ffdc in Eq.(2.3) is nothing but the Susskind lattice fermion and n/ is the number of the 
species doublers (n / = 2) . Restoring the local gauge invariance by Ai and imposing local constraint 
flf x = 1 by introducing a gauge field A , an effective theory for low-energy excitations is given as 
follows in the continuum, 



where the notations are standard and the ellipsis in Eq.(2.5) represents higher-derivative terms. We 
have set K = 1 and also the "speed of light" (the speed of V'-quanta) v s = 1. The original AF 
Heisenberg model is the half-filled state of f x because the number of f x is unity at each site. In the 
effective field theory of the Dirac fermions, this half-filled state just corresponds to the Dirac sea in 
which all the states with negative energy are filled. 

Besides the original spin SU(2) symmetry, the continuum action S^mon has the relativistic in- 
variance as well as the flavor symmetry. It is interesting to see how the original spin operators 
are expressed in terms of the Dirac fermions. From the discussion in Refs.[9, 3], the Neel vector 
N x = {-) X1+X2 S X is related with V as 



where a is the Pauli spin matrix acting on the spin index. 

If nonperturbative effects of the gauge field are negligible and the system (2.5) is in a de- 
confinement phase of the gauge dynamics, the system can be studied, e.g., by the 1/n/ expansion. 
Composite operators of ip may acquire an anomalous dimension, but the massless Dirac fermions 
survive as low-energy excitations at long distances. However if the system is in a confinement phase, 
the structure of the groundstate and low-energy excitations are drastically influenced by A^ . For 2D 
case (and probably also for 3D strong-coupling case), it is expected that for sufficiently large n/ the 
gauge system is in a deconfinement phase even for the strong-coupling limit without the Maxwell 
term, whereas it is in a confinement phase for small n/. That is, there exists a critical number of the 
massless field separating confinement and deconfinement phases. Very recently we estimated this 
critical number for massless two-component bosonic systems by investigating multi-flavor CP 1 
models in (2 + 1)D and obtained n c j = 14, which is rather large compared to the spinon number in 
the typical ASL's. Then it is interesting to see how the groundstate and the low-energy excitations of 
the nodal spinons are changed if the nonperturbative effects like instanton condensation occurs and 
the gauge dynamics is in a confinement phase. We address this problem in the following sections. 




(2.5) 



X1+X2 



(2.6) 



•x 
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3 Lattice Gauge Theory of Spinons 



In this section, we shall study the gauge system (2.5) assuming that the gauge dynamics is in the 
confinement phase, i.e., n/ < nfj. We consider the system at finite-temperature (T) and with finite 
density of spinons. Motivation to consider the finite-density system comes from not only the hole- 
doped high-T c cuprates but also from the physics of the heavy- fermion materials. The Hamiltonian 
of a Kondo-Heiscnbcrg model for the heavy fermions is given as, 

Hkh = ^2 e kc\ a c ka + Jk^Sx- c x ac x + H AF , (3.1) 

k x 

where Cka represents the conduction electrons with the energy efc, and Jk is the Kondo coupling. 
The Kondo coupling prefers hybridization of the conduction electron c x and the localized spin f x , 
whereas the AF coupling tends to make the AF order for S x . In the real experiments, external 
conditions like an external pressure enhances the itinerant properties of the (localized) electrons. 
This fact implies that the external pressure changes density of states and effectively increases the 
Kondo coupling Jk and as a result the number of electrons f x making a spin-singlet pair with the 
conduction electrons is also increased. (See Fig.l) In the following discussion, the density of spinons 
is controlled by introducing the chemical potential /z s instead of considering the Kondo coupling 
explicitly. 




Figure 1: Hybridization of the localized spin and the conduction electron in the AF background of 
the localized spins. Spin-singlet state of f x (a solid arrow) and conduction electron (a dashed arrow) 
appears as a result of the Kondo coupling Jk- 



In order to study the strong-coupling gauge theory for the gauge system (2.5), we introduce a 
e?-dimcnsional spacetime hypercubic lattice. The lattice sites are denoted again by x, but x does not 
represent individual atoms of the original AF Heisenberg model; rather we are considering "coarse 
grained" lattice model valid at long length scales where the gauge-field fluctuations are very strong. 
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Lattice model is obtained by replacing derivative with lattice difference in the action (2.5), 

Sip = \ [^7m c/ m( x )^+m - 4>x+^l P iUl(x)'il) x , (3.2) 

X,[l 

where hereafter fi = 0, • • • , d and U^x) is the compact U(l) gauge field defined on the links (x, x+^i). 
However as well-known in the lattice gauge theory, the action contains additional low-energy 
modes besides <~ 0, the species doublers. In order to decrease these extra low-energy modes and 
to make the symmetry of the system clear, we make the following change of variables, 

i> x =T(x) Xx , i> x =x x T\x), T{x) = { l0 ) x °---{ ld y*. (3.3) 

Substituting Eq.(3.3) into the action (3.2), we obtain 

S 4> = \Y1 [^(^XxU^Xx+fi - f] ll {x)Xx+^Ul{x)Xx , (3.4) 

where r/ (x) — 1 and i]j(x) = (— ) Xo_l ^-i. From the action (3.4), it is obvious that the Dirac 
spinor indices are diagonal in the X x representation, and then we can decrease the number of the 
low-energy excitations by reducing the number of spinor indices. Hereafter let us put the "flavor 
number" of X x as Nf (i.e., Xx, a — 1>" >Nf). In the continuum, the above system (3.4) contains 
2 [(d+i)/2] . Nf masslcss Dirac fermions. Therefore 2[( d+1 )/ 2 l • N f = 2n f = [spin degrees of freedom] 
x [the number of nodes] . 

As we explained above, we shall consider the system of a finite-density spinons. Then the action 
is given as follows by introducing the chemical potential /x s [10], 

5 x = \ [r^(x)XxU li (x)xx+^ - r^Wxx+pUlWxx , (3.5) 

X,[l 

where r (x) = e^ s and rj (x) = r/j (x) . Because of the particle- hole (electron-positron) symmetry, the 
physical results are unchanged under a transformation fi s — > — fi s . We assume ^ s < hereafter for 
the density of spinons decreases by the hole doping or the Kondo coupling. 

Before starting study of the spinon-gauge system, we have to notice the following things, 

1. As we consider the strong gauge-coupling limit, there are no plaquette (no Maxwell) terms of 
Ufj,(x) in the action. 

2. As Uo(x) is the Lagrange multiplyer field for the local constraint on the number of Xx (fx), 
we have to introduce another matter field representing the "background charge" to make the 
total charge vanishing. In real materials, this degrees of freedom correspond to doped holes or 
spin-singlet pairs made of a conduction electron and S x . 
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The requirement 1 in the above is obvious. For 2, we shall introduce the matter field ip x (a = 
1, ■ • • , Nf) by adding the following terms to the action, 

1 



r'o<PxU {x)(p x+6 - (r' ) 1 <P x+ c t U^{x)(p x , (3.6) 



where r' — e Ms . From S v , it is obvious that ip x carry the same charge with Xx with respect the 
Lagrange multiplyer field Uq{x) and they have no hopping terms and are always staying at site x. 
For the original AF system with "hole doping" , the slave-particle representation is useful, in which a 
hard-core boson operator b x is introduced for the "hole state" (or spin-singlet state) and the physical 
state is restricted to (f£f x + b x b x — l)\phys) = 0. See Fig.l. The above <p x on the coarse grained 
lattice corresponds to b x on the original lattice. 2 

In the later calculation of the free energy, we shall see that the chemical potentials fi s and fj,' s 
always appear in the linear combination fj, s — fj,' s as a result of integration over Uo (x) . This means 
that the summation of densities of \x and ip x is automatically vanishing as required by the above 
physical-state condition. 

From S x (3.5), the gauge system has the following global U(iV/)x U(iV/) symmetry; 

Xx^WdXx, Xx^XxWl W D GU(N f ), (3.7) 



and for W A € U(N f ), 



Xx 




x e even site 
x G odd site 



XxWa, x € even site 
Xx - { (3.8) 
XxW A , x £ odd site. 

It is shown that the "Neel order parameter" ij x aip x — Xx&Xx is connected with ip x ip x = XxXx by 
the second symmetry in Eq.(3.8), e.g., 

XxXx -» Xxe^^e^Xx = iXxVjXx, (3.9) 

where cr^'s are the Pauli spin matrices acting on the spin indices and they are the generators of 
SU(2)cU(7V/). For the spinon-gauge system in the continuum (2.5), it is shown that ^ x atp x = Xx^Xx 
and 4) x *l) x = XxXx has the same scaling dimension to all order of the 1/n/ expansion. In the present 
lattice model, these operators are connected by the explicit symmetry. 

In the following section, we shall study the lattice spinon-gauge model in the strong-coupling 
limit by using recently developed techniques of the lattice gauge theory [8]. 

2 More precisely, the flavor number of ip x sholud be Nf/2 for the hole-doped case. In this sense, we are think of 
the two-channel Kondo-Hciscnbcrg model in which two kinds of spin-singlet pair exist. Qualitative phase structure is 
the same for the Nf and Nf /2 cases. 
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4 Free Energy of the Lattice Spinon-Gauge System 

The partition function Z is given as, 

Z = J[DxD X D^DU DU 3 }e- s -- s -, (4.1) 

where j = 1, • • • , d. As we are considering the finite-temperature system, we impose the (anti)pcriodic 
boundary condition on the field variables. Temperature T is given as T — {aN^)^ 1 where a is the 
lattice spacing and No is the number of the temporal sites. Hereafter we simply set a = 1. Unit of 
energy is Tiv s /a and will be set unity. 

There are no gauge-plaquette terms because we are think of the strong-coupling limit. Then we 
can perform the integration over the spatial gauge field Uj(x) (j = 1, • • • , d); the one- link integral. 
We employ the 1 / (i-expansion as in the most of calculations and obtain, 

j [£)C/ . ]e -IE^te(-)^^(-)^-^(-)^ + ^i(-)^] = cxp \j2^(M x V M {x,y)M y )\ , (4.2) 

where 

1 d 

Mf = x a x xl V M (x,y) = -J2(S x ,y + J + S x>v - j ). (4.3) 

Let us introduce auxiliary fields & x b corresponding to the composite fields Jiff. To this end, we 
use the following identity, 

exp \Y^Tv{M x V M {x,y)M y ) 

x,y 

= f[D^}exp[-Y / {^xVM(x,y)^ y ) + 2Tr(^ x V M (x,y)M y )}]. (4.4) 

x,y 

In the following, we shall calculate the free energy (the effective potential) of $ x . Because of 
the symmetry Eqs.(3.7) and (3.8), we can simply set $" fc = & x 5 ab . (See for example, Eq.(3.9).) 
Furthermore we consider homogeneous configuration like § x = $, because it is the lowest-energy 
configuration. Low-energy excitations of <f> x will be discussed in Sec. 6. Then, each term in Eq.(4.4) 
becomes as, 

£ Tr(* x V M (x, y)%) => £ ^% 2 , (4.5) 

x,y x 

and 



J2^ x V M (x,y)M y ) => -^>J2 M x a = 4*E**Xx- (4.6) 

x.y xx 

From Eqs.(4.5) and (4.6), the partition function can be written as 

Z = J[DU }[DxDx}mDip}e- s '^^- s ^ (4.7) 
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where 



S'[U ,x, [lxxe»U (x) Xx+b ~ \x x+ ^ul{x) Xx + ^$ 2 ^XxXx 



(4.6 



Next we shall integrate over Xx and ip x . As they satisfy the anti-periodic boundary condition in 
the 0-th direction, we introduce the Fourier transformed fields, 

1 No N„ 



~'x(n,x) 



(4.9) 



ra=l v " n=l 

where p„ = 7r(2n— 1)/Nq and similarly for (p^. In order to perform the path-integral over the Fourier 
transformed fields, we take the Polyakov gauge for Uq(x) as follows, 



U (x ,x) = e i +W/ No for all x . 



(4.10) 



After integration over x{n, x) etc, we have 

No/2 



x n=l 



(?)' 



Nf 



sin p n {4>) 



Nf $2 



where 



J3„(0) = p n + — lll a , Pn {4>) = p n + — l^ s . 

iV JV 



(4.11) 



(4.12) 



The product ]J^{ 2 in Eq.(4.11) can be calculated by using the usual techniques of the finite- 
temperature field theory, 



Z = J [D<t>] JJ [cosh(iVo£0 • cos(0(x) - iN ol j,' s ] 



1 1 ~i N f dN f ff,2 

+ - cos(20(x) - iN (n s - £)) + - cosh(A ( Ms - ^))J • e ~ * , (4.13) 
where E — arcsinh^^ 

Finally we have to integrate over <f>{x) in Eq.(4.13) to obtain the free energy F e g — — log Z/N V S , 
where the temperature T = 1/N and V s — J2s- Unfortunately, the explicit form of F e g cannot be 
obtained for general Nf. For each Nf, we integrated over 4>{x) in Eq.(4.13) and then F e g is obtained 
as follows for Nf — 2, 



F eS = ^ - Tlog [I cosh 2 (f^A) + \ cosh 2 (§) + \ 



(4.14) 



For N 



f=*, 



F eS = d$ 2 -Tlog 



1 cosh 4 cosh 2 



16 



+ - cosh' 



(f) 



32 



■ cosh 



T 
I T 



(^W(f) 



3 
27 



(4.15) 
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Similar expression is obtained for other values of Nf. As stated above, the chemical potentials 
always appear in the combination yu s — /j,' s and therefore densities of Xx and ip x are automatically 
related with each other, 

dF cB dF cB 

p * = -^7 = ^r = - p ^ (4 - 16) 

as required by the original local constraint f£f x + b\.b x = 1. Hereafter we set \xt = — Ms + Ms > 0. 

In the following section, we shall study the phase structure of the spinon-gauge system by using 
the free energy obtained in this section. To this end, we shall assume a typical form of F e g instead 
of investigating some specific cases separately. 



5 Phase structure of the spinon-gauge system 

As we show in the previous section, the free energy F e g has a rather complicated form. In this 
section we shall study the phase structure of the spinon-gauge system assuming the following "phe- 
nomenological" free energy, 



F eff = ^ 2 -Tlog 



cosh 



(^) +B cosh(^)+c], (5.1) 



where B and C are parameters. For general Nf, there appear more complicated terms as in Eq.(4.15), 
but the essential feature of the free energy and the phase structure can be described by Eq.(5.1). 
First let us consider the small $ region. For $ <C 1, 



2 



4 8 [ cosh (^) + S + C] T V ' 

As long as the coefficient of the $ 4 -term is positive, F cB exhibits the second-order phase transition 
and the chemical potential at the critical points is given as a function of T, 



T 



N f 



BdNf 



P-t = -r^-arccosh — B — C . (5.3) 



2T 



In the region [i T < the field <& has a nonvanishing expectation value and the symmetry is 
spontaneously broken. In the AF spin model, an AF order exists in this region. 

In particular from Eq.(5.3), the critical T at [It = 0, which corresponds to the pure Hcisenbcrg 
model without hole doping, obtained as 

_ BdNf 

Tn ~ 2(1 + B + Cy (5 - 4) 

which corresponds to the Neel temperature. 

As T is lowered, the coefficient of $ 4 becomes negative and the phase transition changes from a 
second-order to first-order one. We call this point tricritical point. The tricritical temperature T tr j 
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can be easily obtained from Eq.(5.1) as 



N f 

Ttn = -f- (5.5) 



From Eq.(5.3), the chemical potential at the tricritical point is given as, 

pf = I • arccosh((d - l)B - c) . (5.6) 

If b +q +1 < 1, the phase transition is everywhere of first order. 

Next let us consider the T = case; the finite-density quantum phase transition (QPT). In this 
case, the free energy reduces to, 

F eB m = ^-$ 2 - JV> • max{/XT, (5.7) 

At $ = 0, E[<& = 0] = and therefore F e g[& = 0] = -Nfp T . It is obvious that <I> = is a local 
minimum of F e fj [<&]. There is another local minimum at $o, where $o is the solution to the equation 
5F e ff[$]/5$ = for E[<&] > p T . <J> is easily obtained as 



and the local minimum is calculated as 

F cS [*„] = - N f E[$ ] < 0. (5.9) 



As ht increases, the global minimum changes from $ to <f> = 0. The critical chemical potential is 
then given by 

p c T0 = E[* ] - \< (5-10) 
and the expectation value of $ behaves as 

{$o for Mt < P'rn 
(5.11) 
for px > Pro- 

Similarly, the density of the spinons p = p x — —p v is calculated from Eq.(5.7), 

dF eS J for p T < Mto ,, 10 * 

P=-5 = < (5-12) 

^ T ( -JV> for p T > p c T0 . 

The above result Eq.(5.11) and (5.12) indicates that a strong first-order phase transition from 
ordered to disordered phases takes place at T = 0, and an abrupt change of the density of the 
spinons from the empty state to the saturated state also occurs there. For px > Mtc the density of 
the localized electrons participating in AF magnetism are vanishing i.e., all electrons are itinerant. 
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Finally let us show numerical calculations for some typical values of the parameters in F e ff. We 
put d = 3, B = 1, C = and N f = 4 in Eq.(5.1). In this case, Tjv = 3, T tri = 2, /4 ri = 
0.658 • • • , fij, = 0.549 • • • and <f>o = 0.693 ■ • •. In Fig. 2, we show the phase diagram in the \xt — T 
plane. The spinon-gauge system at the strong-coupling limit has the symmetry-breaking phase in the 
low-T and low-density region. In the high-T region T > T tr i, the phase transition is of second-order 
whereas in the low-T region T < T tr ; it is of first-order. At quantum critical point (/ir, T) = (/i™, 0), 
the density of spinons changed drastically from zero to —Nf. The phase transition line show a 
reentrant behavior, but this behavior may be an artifact of the mean-field approximation. See 
discussion in the following section. In Figs. 3 and 4, we show the expectation value of $ as a function 
of [It- For T < T tr i, 3> exhibits a discontinuity at (j,t — Mto- Figs. 5 and 6, we show the density p 
as a function of /zt- 




Figure 2: Phase diagram obtained from the free energy by using MFT like approximation. Low-T 
and low-density region of spinons has a AF long-range order. Symbol • denotes the tricritical point 
at /xt = Ht 1 - Above (Below) Tin, the phase transition is of second (first) order. 



6 Discussion and conclusion 

In this paper we studied the strong-coupling region of the spinon-gauge system on the lattice at 
finite-T and finite density of spinons. By using the MFT like approximation, we found that at low-T 
and low density the system has the AF long-range order (LRO). We expect that this LRO is weak 
magnetism compared with the microscopic local Neel order as explained in Ref. [11] and observed 
by the experiments of the heavy-fermion systems. Low-energy excitations in the ordered phase are 
easily identified as follows by the symmetry of the system Eq.(3.8) and also from the result of the 
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Figure 3: Expectation value $ as a function of the chemical potential fix — — Ms + Ms at 7 1 = 1.5 < 
Tin. The phase transition is of first order. 




Figure 4: Expectation value $ as a function of the chemical potential /it = — Ms + Ms a t T = 2.5 > 
Ttrj. The phase transition is of second order. 
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Figure 5: Density of spinons p(< 0) as a function of the chemical potential fix — — Ms + Ms &t 
T = 1.5 < T tri . 
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Figure 6: Density of spinons p(< 0) as a function of the chemical potential fix — — Ms + Ms a t 
T = 2.5 > Ttrj. 
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one-link integral Eq.(4.2), 

^^(^•exp^^rV^)], (6.1) 

K 

where e(x) = (— ) x °^ ^ x d-i j j^k j g ^ ne generators of U(iV/) and ir K (x) are gapless excitations, the 
spinwaves. Then one can easily expect that an effective field theory for the ordered phase is a XJ(Nf) 
nonlinear tr-model. 

At highcr-T, the effect of fluctuations of ir K (x) is larger and destabilizes the AF order. Then the 
reentrance of the transition line observed in the MFT like approximation may disappear as a result 
of the fluctuations described by ir K (x). 

As discussed rather in detail in Ref. [12], the MFT- like approximation is reliable for small Nf, 
whereas for large Nf an instability of the state with long-range order to a symmetric state appears. 
Therefore the strong-coupling gauge theory and also the MFT-likc approximation employed in the 
present paper is justified only for small Nf cases. 

In this paper we considered the gauge theory of spinons. It is quite interesting to study a 
coupled system of the conduction electrons and spin degrees of freedom and discuss the quantum 
phase transition to heavy-fermion state[ll]. 
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